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Why the nature needs 1/f-noise 


Yu. E. Kuzovlev 

Donetsk Free Statistical Physics Laborator'^ 

Low-frequency 1/f-noise occurs at all levels of the nature organization and became an actual factor 
of nanotechnologies, but in essence it remains misunderstood by its investigators. Here, once again 
it is pointed out that such the state of affairs may be caused by uncritical application of probability 
theory notions to physical random phenomena, first of all the notion of’’independence”. It is shown 
that in the framework of statistical mechanics no medium could provide an inner wandering particle 
with quite certain value of diffusivity and mobility, thereby producing flicker fluctuations of these 
quantities. This is example of realization of universal 1/f-noise origin in many-particle systems: 
dependence of time progress of any particular relaxation or transport process on the whole system’s 
detailed initial microstate. 
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I. INTRODUCTION 

A. Root of the question and popular hypothesis 

Here, as for many years before, the question of 1/f- 
noise grows in urgency, by extending and deepening to¬ 
gether with physical experiments and new technologies 
and concerning almost all in the world, from cosmic phe¬ 
nomena down to molecular biology and nano-electronics. 
Nevertheless, there are no modern reviews of the question 
proportional to its volume and significance. Seemingly, 
this is so because investigators do not find an inspira¬ 
tional ideas and happy thoughts for that. Though, one 
appropriate suggestion, - to be under consideration be¬ 
low, - was made already in CHI , but it had not excited 
a response. Somehow or other, today we see reports on 
more and more inventive and fine measurements of 1/f- 
noise, - for example, in films of metals and alloys @ or 
atomic layers of graphene @ , - but as before with no un¬ 
ambiguous indication in its origin. Citing 0, “... despite 
almost a century of research, 1/f noise remains a con¬ 
troversial phenomenon and numerous debates continue 
about its origin and mechanisms”. 

It can be added that “debates”, in the author’s expe¬ 
rience, not rarely take rather totalitarian forms. Maybe, 
partly by this reason, from the author’s viewpoint, the 
present situation in general very slightly differs from what 
was outlined in Cl and a little later in 0] . We would like 
to compare it with situation in astronomy nearly three 
hundreds years ara before appearance of the celebrated 
I. Newton’s work Q. 

We venture such the comparison not for the sake of wit¬ 
ticism but in view of our intention to demonstrate in the 
present notes that just the Newton’s laws of mechanics 
may be the place where solution of the 1 /f-noise problem 
is hidden. More precisely, 1/f-noise is permanent prop¬ 
erty of systems of many particles moving and interacting 
by these laws (in their classical or quantum formulation 
including fields). In order to recognize it, we have only 
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to follow the Newton’s advice to avoid unnecessary hy¬ 
potheses (“Hypotheses non lingo” 0). 

Well, what hypotheses are thought up by physicists in 
respect to 1/f-noise? Let us decipher it by example of 
electric current noise in a conductor under fixed voltage. 
Presence of 1/f-noise there means that the current has 
no certain value, in the sense that its averaging (smooth¬ 
ing) over time produces an unpredictable result, - that is 
randomly varying from one experiment to brother, - with 
a diversity which practically does not decrease, or even 
increases, when the averaging duration grows (since re¬ 
lated narrowing of frequency band contributing to the di¬ 
versity is almost, or with excess, compensated by growth 
of the noise power spectral density inside that band). So, 
when asking oneself a question about origin of such the 
phenomenon, one first of all assumes that it is in some 
specific fluctuation processes influencing the current, - 
through e.g. number of charge carriers or their mobil¬ 
ity, - while specificity of these processes is in extremely 
wide variety of their time scales (memorw or life, or re¬ 
laxation, or correlation times, etc.) Just this is 

main hypothesis. 


B. Idea of the answer and plan of doing 

Really, this hypothesis is not necessary, since the me¬ 
chanics as it is in no way requires certainty of the cur¬ 
rent and, hence, some special reasons for its uncertainty. 
Indeed, no matter what a concrete mechanism of con¬ 
ductivity may be, if it is indifferent in respect to amount 
of charge early transported through the conductor from 
one side of an outer electric circuit to another and thus 
to past value of time-smoothed current, then later this 
mechanism also will be indifferent in respect to them, 
and on the whole it will be unable to set conditions for 
certainty of the current. Thus, it by itself serves as mech¬ 
anism of 1/f-noise. 

What is for the indifference, it is supported by the 
experiment conditions in themselves which state that 
fluctuations of (time-smoothed) current do not meet a 
back reaction of outer circuit instead passively swallow¬ 
ing them. 

In this reasoning, there is no collections of large char¬ 
acteristic times, instead a single time only is present, 
- usually small in practice, - which indicates ending of 
memory of conductance mechanism. If, for instance, it is 
less than several hours, then transfer of charge carriers, - 
with their collisions, scatterings, reflections, etc., - now, 
at present time interval, is passing indifferently to what 
amount of charge was transported yesterday, even if an 
experimental device was not switched off before going to 
bed. Correspondingly, at frequencies lower than inverse 
day one can find a 1 /f-noise. 

Analogously, if somebody possesses unlimited possibil¬ 
ities of profits and expenses and does not keep count of 
them, then he himself could not know how much his ex¬ 
penditures may be on average over time, and it can be 


expected that they will be distributed in time like I/f- 
noise. 

If, returning to the conductor, we short out it, then 
the current’s I/f-noise disappears along with directed 
current. But irregular charge displacements in oppo¬ 
site directions do continue, at that again indifferently to 
their past amount, and thus to their time-average inten¬ 
sity. The latter therefore is not aimed at a certain value, 
which results in 1/f-fluctuations of intensity (power spec¬ 
tral density) of thermodynamically equilibrium “white” 
(thermal) current noise. They are connected to the I/f- 
noise in non-equilibrium current-carrying conductor by 
means of the “generalized fluctuation-dissipation rela¬ 
tions” minsi. 

If one measures equilibrium thermal noise of poten¬ 
tial difference between sides of opened conductor, - e.g. 
electric junction, - then I/f-fluctuations of intensity of 
this noise can be found too. They say that sum of num¬ 
bers (per unit time) of random charge carrier transitions 
from one side to another and backwards is not tracked 
and regulated by the system, in contrast to residial of 
that numbers [1^. At that, characteristic time constant 
of the system (equivalent RC-circuit) determines upper 
time scale for fluctuations of the residial and lower one for 
fluctuations in the sum (while their upper time scale does 
not exist, since they do not change system’s macrostate). 

The aforesaid can be easy extended, - under non¬ 
principal substitutions of particular terms and meanings, 
- to other manifestations of 1/f-noise in the nature. A lot 
of various examples was exposed in BiiliaMElIIl. 
Our demonstration below will be realized in terms of 
equilibrium “molecular Brownian motion” 0, Es m, 

liiMlll. 

In a maximally simple way we shall show that as¬ 
sumption that Brownian particle obeys a certain diffu- 
sivity (rate, or coefficient, of diffusion) is incompatible 
with exact equations of statistical mechanics, that is with 
mechanical dynamical background of Brownian motion. 
Thus, mechanics inevitably generates I/f-, or “flicker”, 
fluctuations of diffusivity and mobility of the particle. 

Then we shall consider quantitative characteristics of 
this I/f-noise and, finally, present its explanation in 
the language of theory of deterministic chaos in many- 
particle systems. 


II. PHENOMENOLOGY OF BROWNIAN 
MOTION 

A. Formulation of problem 

Let us imagine a small “Brownian particle” in a three- 
dimensional statistically uniform isotropic and thermo¬ 
dynamically equilibrium medium. Very small particle of 
dust or flower pollen, - whose motion in liquid for the first 
time was observed through microscope in (^ . and in 
the beginning of next century theoretically analysed in 
(2^ - 1?^ . - are suitable objects. But it will be better to 
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take in mind some “nano-particle” or even merely sepa¬ 
rate atom or molecule in liquid or gas [2^. In principle, 
we may speak even about free charge carrier or point de¬ 
fect in a solid, but confine ourselves by a particle which 
quite definitely is subject to the classical variant of me¬ 
chanics. 

Let R{t) and V{t) = dR{t)/dt denote vectors of 
centre-of-mass coordinate and velocity of our Brownian 
particle (BP) at given time instant, while R and V their 
possible values. We can think that initially at time t = 0 
BP was placed at definitely known space point. Where 
namely, is of no importance, because of thermodynam¬ 
ical equivalence of any BP’s positions. Therefore it is 
convenient to choose the coordinate origin: R(0) = 0. 
Then later instant current position of BP, R(t ), will be 
coinciding with vector of its total displacement, or path, 
during all previous observation time. 

Now, let us ask ourselves what is BP’s “diffusion law”, 
i.e. probability distribution of BP’s path. Density of this 
distribution will be designated by W {t, R). It can be 
represented by expression 

Wit,R) = {6iR-Rit))), (1) 

where the Dirac delta-function figures, R{t) is thought 
as result of all the previous interaction vetween BP and 
the medium, and the angle brackets designate averaging 
over the equilibrium (Gibbs [l^) statistical ensemble of 
initial states of the medium and initial values of BP’s 
velocity. 

Undoubtedly, a plot (relief) of W {t, R) as function of 
R looks like a “bell” extending and lowering with time. 
We are interested in what shapes of this bell may be 
formed in reality. 


Obviously, ([2]) is “continuity equation” for the proba¬ 
bility density W {t, R ), and the “field of velocity of prob¬ 
ability flow”, V(t, R ), containes important information 
about solutions to this equation. Therefore, first of all 
let us consider possible constrution of the vector-function 
V{t,R). 


C. Conditional average velocity of Brownian 
particle 

We shall keep in mind that the duration t of our ob¬ 
servations of BP is much longer than characteristic re¬ 
laxation time r of (fluctuations of) BP’s velocity. 

Then, firstly, apply heuristic reasonings as follow. On 
one hand, by the condition R(t) = R , average value 
of BP’s velocity in the past, at time of its proceeding 
observation, appears equal to R/t. On the other hand, 
as far as BP makes a random walk and t ^ t , the same 
condition R{t) = R tells us almost nothing about BP’s 
velocity in the future, so its average value in equal next 
time interval can be expected to be zero. Hence, since the 
average under question, V (t, R ), relates to the present 
time instant “in the middle between past and future”, it 
seems likely that it is equal to half-sum of the mentioned 
quantities: 

V{t,R) = ^^. (3) 

We can confirm this conclusion in a more formally rig¬ 
orous way, basing on the main distinctive statistical prop¬ 
erty of Brownian motion : 


B. Conditional averaging and continuity equation 

In fact, m is mere identity, but its time differentiation 
immediately brings us a food for thought. From it we 
have 


^EM = _v.(u(t)d(i?-i?(t))) 

(• will denote scalar product of vectors). By attracting 
mathematical tools of the probability theory [s^, this 
equality can be rewritten as 


dW{t,R) 

m 


V ■V{t,R)W{t,R), 


( 2 ) 


where V{t,R) = (U(t))_R is conditional average value 
of BP’s instant velocity determined under condition that 
its current position, and thus its previous path, is known 
(measured) to be equal to R{t) = R. Generally the 
operation of conditional averaging ■ ■)r is defined by 
formula 


{...)R = {...6m)-R))/{6{R{t)-R)). 


{R^{t)) = J R^W{t,R)dR = 6Dt (4) 

at t T, i.e. ensemble average of squared BP’s dis¬ 
placement grows proportionally to observation time. It 
is sufficient to notice that the continuity equation implies 

^ J R^WdR = 2 J R-VWdR 

and that this requirement is naturally satisfied together 
with (jH) (with taking account of parallelity U || i?) when 
equality ([3]) is valid. 

By way, notice that BP’s diffusivity, or diffusion coeffi¬ 
cient, D and relaxation time r always can be connected 
via relation 


D = V^^t = 


T 

M 


where T is temperature of the medium, M is mass of 
BP, and Vq = ^/T/M its chracteristic thermal velocity. 
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D. General form of probabilistic law of diffusion 
and uncertainty of its coefficient 


After inserting function ([3|) into (ED, one comes to par¬ 
tial differential equation 

dW 

2t—=-3W - R-VW, (5) 

which clearly indicates scale-invariant character of its so¬ 
lutions. We are interested in isotropic (spherically sym¬ 
metric) solutions looking as 

W{t, R) = (2Dt)-3/2 (6) 

with some dimensionless function 'I'(z) of dimension¬ 
less argument z = R‘^/2Dt. In our context it, 
since representing probability density, anyway should 
be non-negative and satisfying normalization condition 
/ WdR = 1 in company with equality o, which surely 
can be done. Then ([61 is most general law of diffusional 
random walk, when typical BP’s displacements are pro¬ 
portional to square root of observation time: R^{t) oc t. 

In particular, taking 'I'(z) = (27r)“^/^ exp(—z/2), one 
obtains the commonly known Gaussian diffusion law, 

W = Wait, R) = exp {-R'^/ADt). (7) 


The corresponding walk is much pleasant for users since 
in rough enough, in comparison with t , time scale its 
successive increments are mutually statistically indepen¬ 
dent. Owing to this, the only parameter of such random 
walk, - its diffusion coefficient, or diffusivity, D - can 
be unambiguously determined from observations of any 
its particular realization, by means of long enough time 
averaging. 

However, similar observations and time averaging of 
non-Gaussian random walk, obeying some of general type 
distributions ©, will produce every time different values 
of diffusivity [ll-lj, 0, [l3l • Indeed, their coincidence, that 
is convergence of all results of time averaging to one and 
the same value, would be impossible without statistical 
independence of increments (at least mutually far time- 
distanced ones) which in turn would mean, in accordance 
with respective limit theorem of the probability theory 
(the “law of large numbers”), that at t ^ t probability 
distribution of total path tends to the Gaussian (“nor¬ 
mal”) one [sTj) . 

This becomes quite obvious if distribution (|6D is rep¬ 
resented by linear combination of Gaussian “bells”: 

/A \ dA 

Such expansions naturally arise in the microscopic theory 
[iM flil - llq . Correspondingly, in place of 'k(z) in © we 
can write 




exp(-z/2C) 
Jo ( 203/2 


uiCOdC- 


( 8 ) 


Function t/(C,0 here plays role of probability distribu¬ 
tion of C = A/U, i.e. random diffusivity of BP A ex¬ 
pressed in units of its mean diffusivity D. The latter 
is formally defined by equality (|3D, while practically one 
may try to determine it with the help of averaging over 
many experiments or many copies of BP. 

The additional argument ^ in this expansion, - if intro¬ 
duced e.g. as ^ = r/t under convention 4'(z, 0) = 'I'(z), 
- allows to take into account violation of ideal scale in¬ 
variance of random walk at ^ 0. First of all, far 

on “tails” of diffusion law, where R^ > is 

z > 1/^. There rate of diffusion achieves values of rate 
of free flight, A ^ Vjt = D/^ . 

Clearly, a correction of tails of diffusion law may 
strongly influence its higher-order statistical moments 
and cumulants, even in spite of ^ <C 1. At that, nev¬ 
ertheless, shaping of IF(t,i?)’s bell in the main stays 
almost unchanged. Accordingly, a change of the function 
V{t,R )required by equation © and condition Q, - is 
as small as ^ is, so that the expression ([3D remains right. 

Notice that the very possibility of long-term viola¬ 
tion of scale invariance automatically presumes non- 
Gaussianity of diffusion law, since Gaussian statistics 
merely gives no place for it (since it would contradict 
the condition (|3D). Already this fact gives evidence that 
Gaussian law is not completely adequate reflection of re¬ 
ality, although in mind of scientists it is firmly associated 
with diffusion of physical particles. At the same time nei¬ 
ther general reasonings leading to ED and ED nor equa¬ 
tion ED by itself in no way dictate the special Gaussian 
choice. Therefore, it is desirable to discuss other possi¬ 
bilities and search for criteria of choice among them in 
the framework of statistical mechanics. 


III. MICROSCOPIC APPROACH 
A. Newton equation and Liouville equation 

Further, let us go from kinematics of Brownian mo¬ 
tion to its dynamics and directly consider BP’s interac¬ 
tion with medium using methods of statistical mechanics. 
With this purpose we can take for our system quite usual 
simple Hamiltonian 

+ (9) 

where P = MV is BP’s momentum, F is full set of 
(canonical) variables of the medium, $(i?, F) is energy 
of BP-medium interaction, and Hth{^) is Hamiltonian 
of medium in itself (or, in other words, that of “thermo¬ 
stat”). If BP possesses internal degrees of freedom, then 
their variables will be thought included into the set P, 
thus being formally treated as a constituent of medium. 

Let V = V{t, R,P,T) designate density of full proba¬ 
bility distribution of microstates of our system. Its evolu¬ 
tion is described by the formally exact Liouville equation 
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[sj. Here we can display it partly, writing out its 
terms only directly concerning BP: 


and the double over-line marking tensor (matrix) of con¬ 
ditional quadratic cumulants (second-order cumulants) 
of velocity: 


— =-V-WV - F{R,T)-WpV + ... . (10) 

Here F{R,T) = — V$(i?, P) is force acting onto BP be¬ 
cause of its interaction with medium, and the dots sur¬ 
rogate terms with P derivatives. 

Considering probability distribution of displacement 
(coordinate) of BP, 


W{t,R) 


J j V{t,R,P,r)dTdP, 


from equation cni) after its integration over P and P 
one comes, of course, to the continuity equation ([2). The 
same integration after multiplying ([T0|) by V produces 
additional equation 


g _ _ _ 

— VW = -S/-VoVW + M-^FW. (11) 

Here and below the symbol o denotes tensor product of 
vectors, while the over-line means, as before, conditional 
averages under given R{t) = R. Namely, in the first 
term on the left 


VoV{t,R) 


{V{t) o V{t))R = 


ffVovvdrdP 

w 


and in second term there 


F(t,R) = (F(R(t),r(t))}H 


ff F(R,r)vdrdP 

w 


Equation CD describes momentum exchange between 
BP and medium. In essence, - as it can be easy verified, 
- this is merely the Newton equation M dV/dt = F after 
its conditional averaging: 


{M dV{t)/dt - F{R{t),r{t)))R = 0 . 


We shall transform it into relation between functions 
F{t,R) and W{t,R) which is able to help selection of 
acceptable diffusion laws without more deepening into 
the Liouville equation. 


B. Equation of friction of Brownian particle 


Replacing derivative dW/dt in the equation (fTTI) with 
right-hand side of Q, after simple manipulations one 
comes to equivalent exact equation 

dV V • VWW F 

W “M’ ^ ^ 

with “material derivative” of BP’s average velocity. 


dt 




+ (y-v)i/. 


V oV = V oV -V oV. 


Next, at first let us consider the latter object. 

Since we are speaking about thermodynamically equi¬ 
librium Brownian motion, we can state that the condi¬ 
tional cumulants’matrix VoV{t,R) at f ^ T coincides 
with matrix of unconditional equilibrium quadratic sta¬ 
tistical moments of velocity, {V(t) o V{t )), that is re¬ 
duces to scalar number Vq = T/M regardless of R. 
Indeed, if t ^ t , then at any R the condition R{t) = R 
fixes BP’s position occurred after many random steps 
and cycles of momentum and energy exchange between 
BP and medium under detail balance in this process. 
Therefore, the value (variance) of corresponding thermal 
randomness of BP’s velocity is not affected by this condi¬ 
tion (otherwise, thermal kinetic energy of BP, on average 
equal to MV -1^/2, would be dependent on where BP is 
found). 

The said can be confirmed by direct calculation of the 
matrix VoV{t,R) in case of Gaussian random walk sub¬ 
ject to distribution 0, which yields 

VW{t, R) = (1 _ ^/2) ^ 1/^2 (13) 

at^ = r/t—>-0. This result is valid also in case of 
non-Gaussian walk obeying dH) and ([8|), since its differ¬ 
ence from Gaussian one is contained in its higher-order 
cumulants. 

Then, let us compare two terms on left side of m- 
For the first of them insertion of expression gives 


^ ^ T R 

H ~ “ 4^2 ~~2M ^ ■ 

For the second term after insertion of m and (l6|) we 
have 

/ A T 2d\n^iz,^) R T R 

~2)m di ^ 


with same shortened notation z = R^/2Dt as before. 
Right-hand expression here corresponds to the Gaussian 
diffusion law, for which dln^'/dz = —1/2, but by order 
of magnitude it is true in general case too, at least at z <C 
1/^. It shows that the first term, being approximately 
2</t times smaller than the second, is negligibly small in 
the limit ^ ^ 0 . 

Hence, consideration of long enough time intervals 
leads us from m to shortened relation 



2dln^'(z,C) 

dz 



(14) 


It resembles an equation of viscous friction, with R/2t = 
V in the role of velocity of a body moving through fluid 
and the round brackets in the role friction coefficient. 
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One more simplification can be obtained by neglecting, 
under mentioned limit, violation of scale invariance and 
treating as a function of single argument '^{z). 

In the next paragraphs we firstly proceed just so. 

But before that let us once again glance at the van¬ 
ishing first left term of lfT^ . If writing its contribution to 
the mean force as 

,, dy „ My^ 

M — = - V —— , 
dt 2 

one can say that this is force of reaction of the medium 
_ 2 

to addition MV 12 to energy of BP, and system as the 
whole, introduced by the very measurement of BP’s path, 
and therefore this force is not sensible to shape of diffu¬ 
sion law and thus to concrete peculiarities of medium. 
There is evident analogy with perturbing effects of mea¬ 
surements in quantum mechanics. 

In opposite, the remaining, in the large-time limit, part 
of the force, (O, is determined solely by shape of proba¬ 
bility distribution of equilibrium Brownian displacement 
(“diffusion law”). Consequently, this force characterizes 
inherent, - unperturbed by observations, - BP-medium 
interaction. In particular, it shows characteristic levels 
of the interaction forces and energies necessary for real¬ 
ization of one or another concrete diffusion law. Now, 
examine in this respect the Gaussian law m and make 
sure that it is unrealistic. 


C. Paradox of Brownian motion: Gaussian 
statistics for it is beyond strength of its mechanics 


For Gaussian diffusion law, the square bracket in the 
“friction equation” (ITTl) turns to unit, and the equation 
becomes linear: 


F 


T R _ R ^ T 


(15) 


At that, the “friction coefficient” in front of R/2t = V 
connects to the diffusivity via relation similar to the 
widely known “Einstein relation” [2^, [2^. Suchlikeness, 
however, is not a plus but minus of equality (IISl). 

The matter is as follows. Friction force in the true Ein¬ 
stein relation represents medium resistance against di¬ 
rected motion of a particle. When this particle displaces 
by distance R , the corresponding force makes work 


thus producing a heat (recall that z = R^/2Dt). This 
quantity, - like the force itself, - in principle can be ar¬ 
bitrary large under proper initial value of the particle’s 
kinetic energy. 

This is clear. But it is strange thing that equality 03) 
offers the same, also unbounded, characteristic values of 
force and work. Such a picture categorically contradicts 
to sense. 


Really, - repeating the aforesaid, - in our case the force 
what figures in (TTSl) represents medium reaction to par¬ 
ticle’s displacement achieved along random trajectory of 
thermal motion, when initial energy value knowingly is 
only ^ T. At that, the medium creates obstacles to in¬ 
ertial free flight of BP but in no way to its unrestricted 
moving off from beginning of its path. In opposite, the 
moving off proceeds due to medium’s own free will and 
at the expense of its own equilibrium fluctuations. 

Therefore in reality, in contrast to m, the average 
force 0, as a function of passed path R , can not be 
arbitrary large, instead staying always and everywhere 
bounded. This is required by such factual inherent prop¬ 
erty of Brownian motion as translational invariance, that 
is indifference of the system in respect to irretrievable de¬ 
partures of BP anywhere. Moreover, on this ground it is 
reasonable to expect that at large \R\ the returning force 
vanishes at all. 

Thus, we have to conclude that the Gaussian law is 
inadequate to physical origin of Brownian motion. 

Inevitability of this conclusion catches eye when notic¬ 
ing that if equality (HU was true then it would mean 
that medium returns BP to start of its path with force 
proportional to the path, F rx —R, i.e. like ideal spring 
with potential energy zT /2 oc . Prom physical point 
of view this looks absurdly, since any far BP’s going away 
is permitted just because it does not change thermody¬ 
namical state of the system. 

Our conclusion can be denominated as paradoxical, if 
recollecting that Gaussian diffusion law many times is¬ 
sued from pen of theoreticians in various physical con¬ 
texts and occupies important place in idealized world 
of “mathematical physics”. But the paradox resolves in 
very simple way: the Gaussian statistics always had ap¬ 
peared as consequence of clear or implicit hypotheses (or 
postulates) about “independences” of random events or 
values. What is for us, we have managed without such 
hypotheses and thus showed their fallacy in application 
to Brownian motion. 

In the past, we too were not connected with them 
and arrived to the same paradoxical conclusion, in the 
framework of both phenomenological statistical analy¬ 
sis of diffusion and transport processes EH 113 and 
analysis based on the full hierarchy of Bogolyubov-Born- 
Green-Kirkwood-Yvon (BBGKY) equations 0, [13 H) 
as well as on the base of exact “generalized 
fluctuation-dissipation relations” (FDR) or “dynamical 
virial relations” El [il El, [13, and by other meth- 
ods ^ [13 1^ - including that for quantum systems 

In Section IV we shall again touch on “paradox of inde¬ 
pendence” . And now, in next paragraph, consider exam¬ 
ples of physically correct diffusion law as an alternative 
of Gaussian one. 
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D. Thermodynamics of Brownian motion and 
statistics of large deviations 

From the left expression in m it is clear that the 
boundedness of the force F in general implies relation 

\F{t,R)\ ^ (16) 


It must be underlined, besides, that boundedness of 
the force declared by (ITBl) also logically implies vanishing 
of the force at infinity (excluding border case only when 
Fmax{t) = |F(t,oo)|), so that the medium’s “spring” re¬ 
sists to small “stretching” only and always loses elasticity 
at large stretching. 

The appropriate example of diffusion law satisfying 
(HZl, that is possessing power-law tails, is presented by 


whose right-hand part can be easy guessed for reasons 
of dimensionality. Of course, the symbol ~ here hides 
some dimensionless coefficient which reflects particular 
distinctions of the system and construction of the func¬ 
tion - In'I'(z). 

Comparison between dm) and m shows that in the 
region of “tails” of diffusion law, at z>l, the Gaussian 
law requires to exaggerate the real force value at least 
^ -y/i times, thus thoroughly falsely describing (strongly 
underestimating) probabilities of large displacements of 
BP with z ^ 1. 

In reality, according to (ITSl) . function — In'I'(z) grows 
not faster than oc ^/z , so that (— In ^'(x))/ \fz < oo , and 
consequently decrease of W{t,R) at large |i?| —>■ oo is 
always sub-exponential (anyway, not faster than simple 
exponential, not speaking about “Gaussian”). 

The difference of reality from “Gaussian ideal” be¬ 
comes aggravated when not the force itself only is 
bounded but also value of characteristic energy (work) 
conjugated with this force: 

Aiz) = \R\\F\/2 < Amax = Aioo) - T (17) 


(with the same remark about ^). It is natural expecta¬ 
tion in view of that at result of any walk (any path R) 
the medium takes from BP not more energy than BP had 
been able to take from medium before. 

Boundedness of A{z) implies that of the force, more¬ 
over, implies that under increase of |i?| the force passes 
through a maximum and then decreases down to zero, 
approximately as |F| Ri 2,Amax/\R\ oc T/\R\. This 
asymptotic again is prompted already by dimensionality 
of quantities we give to disposal of statistical thermody¬ 
namics. 

As the consequence, following equality (ITTl) . the tails 
of diffusion law and thus probabilities of large deviations 
(z 1 ) from typical behavior (z ^ 1 ) decrease un¬ 
der growth of |i?| even much slower than in mere sub¬ 
exponential fashion generally dictated by inequality dT51) . 
Now they decrease in a power-law fashion: 

'k(z) oc (z —>■ oo) . 


It is seen after scalar multiplication of m by R , then 
solving so obtained differential equation, which yields 


^'(z) = ^'( 0 ) exp 


A{z) 

Tz 


dz 


«'(z) 


( 3/2 + 77 )! 
(27177)^/2 77! 



(18) 


with free parameter 77 > 0 (factorial x! is standard syn¬ 
onym of the gamma-function r(x + 1)). At that, ob¬ 
viously, Amax = (5/2 + rj) T. The condition 77 > 0 is 
necessary for finiteness of the mean diffusivity in (|T|. 

Such distribution, with 77 = 1, for the first time 
was obtained in [l4| from consideration of Brownian 
motion (“self-diffusion” 0 ) of test, or “marked”, atom 
of a gas. Similar distribution was found for molecular 
Brownian motion in a liquid [313). Though, strictly 
speaking, this is approximation of formally more exact 
but more complicated expressions taking into account, 
among other factors, violation of the scale invariance. 

Formula (1181) turned out to be a reasonable approxi¬ 
mations also for BP whose mass M differs from mass m 
of medium (gas) atoms. At that, various mathematical 
approaches [313 [31 to the BBGKY equations lead to 
identical estimate of the parameter 77 as a function of 
mass ratio, rj = M/m . 

Hence, investigation of complete (infinitely-many- 
dimensional) Liouville equation qualitatively justifiers re¬ 
sults of our semi-heuristic analysis of initial terms of this 
equation. 

We may further lower formal rigor and try to visually 
“by fingers” interpret mathematical connections between 
statistics of Brownian motion and its microscopic mech¬ 
anism. For instance, namely, let H be internal pressure 
of the medium (gas) and the quantity Amax = A(oo) 
be identified with 3T/2 + HH, where ft is gas volume 
forced out by far walking BP from its vicinity, and HH is 
related forcing work. In essence ft represents deficiency 
of BP’s collisions with gas atoms facilitating its far go¬ 
ing away. At that, since position of center of mass of 
the system stays fixed, an effective decrease of gas mass 
near BP, mfln , - with n being mean concentration of 
gas atoms, - just compensate local mass excess M + m 
accompanying current BP-atom collision. From here we 
have ft = (M/m + l)/n and, taking li/n = T for not 
too dense gas, Amax = (5/2 -\- M/m)T. 

Such reasonings, of course, by themselves are rather 
unsafe, but they can be supported by exact results on 
pair man y-p article non-equilibrium statistical correla¬ 
tions [3 13, [3 HSl- particular, that is a theorem 
stating that short-range character of one-time spatial pair 
BP-atom correlation (boundedness of “correlation vol¬ 
ume” ft ) implies long-range (“long-living”) behavior of 
many-time self-correlations in BP’s motion and thus in- 


and finally applying inequality m- 







validity of Gaussian diffusion law (and, reciprocally, va¬ 
lidity of the latter reqmres non-locality of BP-gas corre¬ 
lations in space) [H, . 


E. Uncertainty and flicker fluctnations of diffusivity 

The expansion (|5]) of non-Gaussian diffusion law m 
over Gaussian ones yields for related probability distri¬ 
bution of C = A/D expression 


= (?) “>’(“?)■ 

According to the FDR P, 0, [13 this distribution trans¬ 
mits onto BP’s mobility (at least “low-field” one) and 
therefore can be observed in measurements of diversity 
of “time-of-flight” (time of drift) values of Brownian par¬ 
ticles under influence of external force (for example, in¬ 
jected electrons or holes in semiconductors) fisjl - 

Effects of non-Gaussian statistics were observed also 
directly in equilibrium, by measuring “fourth cumulants” 
(irreducible fourth-order correlations) of electric current 
or voltage noise P. At that, low-frequency, at frequen¬ 
cies / ^ l/t, fluctuations of power spectral density of 
thermal white noise were under investigation, i.e. in 
essence uncertainty and fluctuations of rate of charge 
transfer and rates (coefficients) of diffusion of charge car¬ 
riers. 

In our example (fT8|) - (fT^ . for 77 > 1 , it is not hard to 
obtain 


(i?2)2 V(A )2 ) r]-l 

with (A) = D and {R^) = 6Dt . This formula shows 
not only degree of uncertainty of diffusion rate but also 
defects of approximation of pure scale invariance: di¬ 
vergence of variance of diffusion rate at 77 < 1 and 
pure independence of the variance on duration of ob¬ 
servations. The latter makes the rate fluctuations effec¬ 
tively quasi-static, with spectrum (power spectral den¬ 
sity) S'_d(/) oc S{f) concentrated at zero frequency. 
This is usual result of a simplest (though non-trivial) ap¬ 
proach to 1 /f-noise from microscopic theory [ij, [ 2 l| . 

Undoubtedly, in a more precise theory beyond ideal 
scale invariance [illiS the slow tails of diffusion law 
are somehow “cut off”, at least at ^ (A > 

D/^ ^ V^t in dH])), so that the A’s variance, as well 
as all the higher-order statistical moments of R and A , 
definitely are finite and hardly exceed values correspond¬ 
ing to free BP’s flight: {{R^)^) < {2k + 1)!! 
and (A*’) < {2k + 1)!! {VQt/2)^ . What is for the delta- 
function S{f ), it in a definite way “spreads”, with keep¬ 
ing dimensionality and singularity at zero, into ^ 1 //, 
where ^ replaces some function of In (r/). 

The first of these corrections is easy describable by re¬ 
placing U{(,0) , - for instance, in (ITOl) . - by approximate 


expression U{(,^) ~ U{(,0)E{(^), in which S(0) = 1 
and S(-) in sufficiently fast way tends to zero at infinity. 
Then instead of (IT^ one obtains 'I'( 2 ,^) ft! 'k(z) Q{z ^), 
where scale-invariant factor 'I'( 2 ;) is the same as before, 
- for instance, in (IT^ . - and also 0 ( 0 ) = 1 and 0 (-) fast 
decreases to zero at infinity, thus cutting off the 47(z) ’s 
tail. As the result, the quadratic cumulant of A becomes 
finite even at ry < 1. At once it acquires time depen¬ 
dence, so that the fourth cumulant of BP’s displacement 
increases with time oc , if 77 < 1 , and oc In (t/r) 
at 77 = 1. Gorrespondingly, the quasi-static spectrum 
oc S{f) transforms to “flicker” spectrum. 


SD{f) 


£>2 


1 

Tf. 


1-17 


( 20 ) 


at r/ <C 1 , in particular, to oc 1 //, when 77 = 1 . 

However, at 77 > 1 such correction is insufficient for 
full “spreading” of frequency delta-function, which says 
that scale invariance violation in this case has some an¬ 
other or more complex character. A notion of how else 
it ma y lo ok we can obtain, for instance, if consider [ij- 
0 , [13,113 diffusion law possessing property of infinite di¬ 
visibility in the sense of the probability theory na but 
asymptotically only for ^ = r/t —>■ 0 , since no real 
transport process can be physically divided into infinitely 
small independent pieces. The respective kernel in the 
expansion ([5]) simplistically is describable by formula 


U{C,0 


a{^) exp+ [-(C- Co(6)/c] 

[C-Co(0 + «(e)]^ 


( 21 ) 


where a( 5 ) = l/ln(l/^) = [ln(t/r)]“^, function 

exp_|_(a;) =exp(a:) at cc > 0 and exp_|_(a:) = 0 at 
a: < 0 , function Co (0 is determined by condition (| 1 ]), i.e. 
J(U{(,^)dC = 1 , while c = Tq/Dto with tq and tq 
being minimal space and time scales down to which the 
“infinite divisibility” of random walk is physically mean¬ 
ingful presumes for simplicity that the constant c 

is not too small, c >> a(^)). As it is seen from here, 
at ^ ? 0 expression (1^ turns to U{(,0) = 6{({ — 1). 

Thus, the scale-invariant “seed” of such the diffusion law 
is purely Gaussian, which motivates to name it “quasi- 
Gaussian” In it was considered in detail, includ- 
in^ its generalizations and comparison with experiments 

For tails of the quasi-Gaussian law at z 1 from (ED) 
and ([ 8 |) one can find 


'i'(o,e) 



that is tails satisfy the boundedness requirement m, 
although he on boundary of set of diffusion laws permit¬ 
ted by (HE). And for spectrum of flicker fluctuations of 
diffusion rate, - or, generally, rate of a transport process, 
- the kernel ( 1 ^ yields 
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SdU) 


D^c 

~T 


In 



( 22 ) 


where 7 = —2 . 

Spectrum (1221) results mainly from difference between 
most probable and average (ensemble-averaged) values of 
the rate, namely, Co('?) and 1 in ( 1211 ) in relative units. 
More precisely, (l2^ reflects logarithmically slow decay 
of this difference with observation time: 1 — Co (0 ~ 
a (0 ln(c/a( 0 ). 

The kernel (HU) by its shape is quite similar to m 
(both consist of more or less sharp “wall” on the left and 
comparatively gentle slope on the right), but analogous 
difference in (HU) is fixed. Imparting a time dependence 
to it may be one more, parallel, scenario of spreading of 
spectrum oc 6{f) under improved analytical approxima¬ 
tions of solutions to the BBGKY equations. From our 
point of view, this is practically important problem of of 
statistical mechanics. 

Though, even presently available approximations of 
microscopic theory are able to realistic quantitative es¬ 
timates of 1/f-noise amplitude. Estimates obtained in 
0 and in [ij, [3 different approximations (dH]) with 
7=1 and ( 1 ^ with 7 = 1 or ( 1 ^ with 7 = 0 , re¬ 
spectively), although differing one from another by factor 
ln[l/(T/)], nevertheless, both are in satisfactory agree¬ 
ment with experimental data on liquids and gases [l|, IlSl j , 
with taking into account diversity of these data. 

On the other hand, the scheme of quasi-Gaussian ran¬ 
dom walk rather well predicts or explains level of electric 
1 /f-noise in various systems [IHli [3, El- Since trans¬ 
ported physical quantity there is charge instead of mass 
(in view of smallness of mass of usual charge carriers), 
and interactions of walking charges with medium is es¬ 
sentially long-range, it is not surprising that transport 
statistics there is non-Gaussian in essentially other man¬ 
ner than in case of molecular Brownian motion. Analysis 
of relation of this statistics to a quantum many-particle 
Liouville equation or equivalent “quantum BBGKY hi¬ 
erarchy” , - for e.g. standard electron-phonon Hamiltoni¬ 
ans, - also is actually important problem M- 

At today’s stage of development of statistical mechan¬ 
ics it is useful to state that unprejudiced treatment of this 
science inevitably discovers flicker fluctuations of rates 
of transport processes, even diffusivity (rate of random 
walk) of particle in ideal gas [ 3 , E) ES El and, more¬ 
over, even in the formal Boltzmann-Grad limit (under 
vanishingly small gas parameter) [30l| . 

This fact excellently highlights inconsistency of at¬ 
tempts to reduce 1 /f-noise and related long-living sta¬ 
tistical correlations and dependences to some very long 
memory or relaxation times. And thus it highlights in¬ 
consistency of the underlying opinion that any statistical 
correlations between random phenomena gives up some 
literal or at least indirect physical correlations between 
them. 


In the next Section by means of elementary logics only 
we shall show that in reality in many-particle systems 
just physical disconnectedness of inter-particle collisions 
leads to uncertainty and 1 /f-noise of relative frequency of 
collisions and rate of wandering of each particle. Thus we 
from a new viewpoint shall justify both the general logics 
of Introduction and the following elementary mathemat¬ 
ical analysis of molecular random walk. 


IV. MYTHS AND REALITY OF RANDOM 
WALKS 

A. Gaussian probability law and two meanings of 
independence of random events 

First, recall why the Gaussian law have appeared and 
appears in various theoretical models. This is because 
it naturally comes from assumption of statistical inde¬ 
pendence of BP’ displacements (increments of random 
walk) at non-intersecting time intervals. And, most im¬ 
portantly, because physicists have gotten accustomed to 
identify statistical independence of random events in the 
sense of the probability theory with their independence 
in the sense of their non-influencing one on another. 

Both these circumstances have more than three hun¬ 
dred years history. A history of the Gaussian law had 
taken beginning from the celebrated “law of large num¬ 
bers” [3| discovered by J. Bernoulli who investigated 
statistics of sequences of observations on vicissitudes of 
life or, for instance, coin tossing or playing dice, under 
assumption that unpredictable outcomes of successive 
“random trials” are mutually independent. To be more 
precise, that their probabilities are independent, that is 
joint probability of several random events decomposes 
(factorizes) into product of their individual probabilities. 

Exactly in such the way the (statistical) independence 
is introduced in modern probability theory [33j. But 
there it is nothing but formal mathematical definition, 
and therefore, - as A. Kolmogorov warned in [3|, - de¬ 
duction of this probability property from seeming inde¬ 
pendence of physical phenomena as such is possible only 
as a hypothesis to be verified by experiments. 

In other words, any evidences of independence of phys¬ 
ical random events at every concrete their realization, 
- in the sense, for instance, of absence of cause-and- 
consequence connections between them, - as such can not 
be sufficient ground for declaring statistical independence 
of these events in a set (statistical ensemble) of realiza¬ 
tions (observations). 

Logically inverting this thesis, we obtain that even 
when statistical experiments reveal statistical depen¬ 
dence in an ensemble of realizations of random events, 
this observation does not necessarily mean existence of 
some real interaction of the events. Just such situations 
do occur when one meets 1 /f-noise. 

Hence, identifying of the two meanings of “indepen¬ 
dence” is nothing but fallacy. Unfortunately, it tradition- 
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ally governs relations of physicists to randomness, even 
despite its careful disclosure, - from viewpoint of funda¬ 
mental statistical mechanics, - by N. Krylov more than 
sixty years ago (^ . 


B. Collisions, chaos and noise in system of hard 
balls 

Mathematicians know N. Krylov as one of pioneers 
of modern theory of dynamical chaos. According to it 
[13, [13, for instance, motion of iV > 3 elastic hard balls 
in a box or disks on torus, obeying deterministic laws 
of mechanics, is indistinguishable from a random process 
[13[13- For us here, it is important that statistical char¬ 
acteristics of this process are crucially sensible to ratio 
of its observation time and total number of balls partici¬ 
pating in it. 

More precisely, let us consider role of parameter t/rN , 
where r is mean free path time of a given ball and thus 
characteristic time of relaxation of its velocity because of 
collisions with other balls [3[Hi- At 

t/r > N, 

clearly, number of quantities describing trajectory 

of any particular ball is much greater than number oc N 
of quantities establishing initial state of the whole sys¬ 
tem, so that each particular trajectory contains one and 
the same exhaustive information on the system. More¬ 
over, this information is contained even in any small part 
of the particular trajectory with duration ^ Nt <C t . 

Due to this circumstance, fluctuations in numbers of 
collisions of given ball from any time sub-interval ^ Nt 
to next one behave like statistically independent ran¬ 
dom values, or “white noise” (which is well understand¬ 
able: presence of some relationship or correlation be¬ 
tween them would be recognition of some system’s initial 
state specificity yet non-realized on shorter time inter¬ 
vals, in contradiction to the condition t ^ Nt). Cor¬ 
respondingly, relative frequency of the ball’s collisions 
time-averaged over whole observation time is almost non- 
random, that is one and the same for all balls and all ini¬ 
tial conditions (at fixed full system’s energy, of course), 
while statistics of fluctuations in number and rate of col¬ 
lisions (of given ball) at intervals t ^ Nt obeys the law 
of large numbers, i.e. is asymptotically Gaussian. 


C. Paradox of independence 

Such the picture of chaos of collisions like usual noise 
is quite pleasant for physicists. But we should not forget 
that it had required the condition t/T N establishing 
rigid (detetrministic) non-local in time and space (non¬ 
vanishing at t/r —> oo and spanning all the balls) phys¬ 
ical (cause-and-consequence) inter-dependence between 


collisions. Just at the expense of this dependence, - para¬ 
doxically! - the statistical independence of time-distant 
and space-distant events (collisions) was ensured. 

In other words, interestingly, creation of ideal disor¬ 
der, - with which statistical independence is usually as¬ 
sociated, - needs vigilant underlying control of it and, in 
this sense, global strict order. At this point we invol¬ 
untarily remind how Dront in the B. Zakhoder’s Russian 
translation of the L. Carroll’s “Alice’s adventures in won¬ 
derland” agitated other personages to “fit into strict dis¬ 
order” , or “stand up strictly anyhow”. Along with these 
laughable words, comparisons suggest themselves with 
the mysterious “quantum non-locality” and “entangled 
quantum states”. 

D. Uncertainty and 1/f-noise of relative frequency 
of collisions and rate of diffusion 

However, in the real world it is not simple to mark off 
temporal disorder of random events in so strict way as 
to subordinate it to the law of large numbers. It is not 
simple by those simple reason that real many-particle 
systems are characterized by just opposite ratio of du¬ 
ration of observations (practically achievable in experi¬ 
ments) and number of particles in the system: 

t/r <C N. 

Therefore, the appeal to arbitrary large averaging times, 
so much beloved in mathematical physics, has no factual 
grounds [^ . 

The above inverse inequality is satisfied even for rather 
small volumes of solids and fluids isolated from the rest 
of the world [l^l- All the more, this inequality is true if 
one takes into account physical impossibility of complete 
isolation and hence necessity to include to N particles 
(and generally degrees of freedom) of all huge surround¬ 
ings of a system under interest. And definitely this in¬ 
equality covers objects of the Gibbs statistical mechanics, 
in which number of particles N is not limited, and which 
was under N. Krylov’s critical analysis [s^ . 

Now, number oc tjT of quantities sufficient for de¬ 
scription of observed trajectory of one or another particle 
(ball) all the time stays small as compared with number 
oc of independent causes, i.e. variables of system‘s 
(initial) state, determining the trajectory. 

But averaging over relatively few number of conse¬ 
quences determined by much larger number of causes 
definitely is unable to produce a certain result, since the 
result remains dependent on many unknown free parame¬ 
ters and does not represent all possible variants of course 
of events, all the more can not represent them under some 
certain proportion. Therefore, time averaging of observa¬ 
tions of particle’s motion in any particular experiment (at 
each realization of system’s phase trajectory) inevitably 
brings unpredictably new value of relative frequency of 
the particle’s collisions, all the more, new distribution 
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(histogram) of collisions (or more complex events) in re¬ 
spect to their inner characteristics. In other words, an 
experimenter meets 1/f-noise (see Introduction). 

From here we see that, instead of fabrication of hy¬ 
potheses on relative frequencies or “probabilities” and 
“independences” of events constituting random walks it 
would be better for all that to follow Newton Q and de¬ 
vote ourselves to investigation of equations of (statistical) 
mechanics. 


E. Game of independences and problems of 
statistical mechanics 

Just said is just to what N, Krylov called in his book 
[ 3 ^ clarifying falseness of the widespread prejudices (ci¬ 
tation [^) “... as if a probability law exists regardless 
of theoretical scheme and full experiment” and “... as 
if “obviously independent” phenomena should have inde¬ 
pendent probability distributions”. 

The “full experiment” here means concrete realization 
of system’s phase trajectory considered as a single whole, 
- as an origin of practical observations, - without its ar¬ 
tificial division into “independent” time fragments (thus, 
we in Sections 2 and 3 above have analyzed just a full 
experiment). 

As far as, - at N ^ t/r, - time-smoothed relative fre¬ 
quency, or rate, of a given sort of random phenomena or 
events (collisions of a given particle with others) varies 
from one experiment to another, demonstrating non-self¬ 
averaging, we can not (have no grounds to) introduce 
for such event a separately definite (individual) a priori 
“probability”. This means that all the events occur seem¬ 
ing commonly statistically dependent, since, figuratively 
speaking, all equally are responsible for resulting, each 
time new, rate of their appearance (a posteriori proba¬ 
bility) . This is so in spite of that physically all the events 
are independent, since at N ^ t/r are determined by 
interactions with different groups from total set of N 
particles. Consequently, we come to crash (inapplica¬ 
bility) of the Bernoulli’s law of large numbers based on 
postulate of statistical independence. 

Here, we clearly see another side of “paradox of inde¬ 
pendence” : a true full-value chaos implies infinitely long 
statistical dependences and correlations. 

It is clear also why the molecular Brownian motion, 
being conjugated with such full-valued chaos, does not 
want to go into “Procrustean bed” of Gaussian statistics 
and, all the more, Boltzmann’s kinetics. 


V. CONCLUSION 

Unfortunately, the above underlined popular careless 
ideas of independences and probabilities of random phe¬ 
nomena (once again citing [36|) “... are so much habitual 
that even a person who had agreed with our argumenta¬ 
tion then usually automatically returns to them as soon 


as he faces with a new question. The origin of stable¬ 
ness of these ideas is in that they are based on com¬ 
mon intuitive notion about statistical laws, and therefore 
they would be permissible and advisable if the talk con¬ 
cerned learning of phenomena of empirical reality. How¬ 
ever, such ideas turn out to be quite unsatisfactory as a 
bench-mark for substantiation of probability laws when 
the talk is about connections between statistical laws to 
principles of the micro-mechanics”. 

Fortunately, at present we have understanding of errors 
of replacing micro-mechanics by speculative probabilis¬ 
tic constructions, let beautiful in themselves and likely. 
Besides, as we noted above, there is already an experi¬ 
ence of consecutive investigation of equations of statis¬ 
tical mechanics in application to transport processes. It 
clearly shows that mechanics of systems of very many 
interacting particles, or degrees of freedom, in now way 
prescribes for the interactions to keep definite rates of 
changing system’s micro-state (transition probabilities), 
even when molecular chaos takes form of a macroscopic 
order (let even thermodynamic equilibrium). 

The point is that any realization of “elementary” act 
of interactions in fact is a product of full (initial) micro¬ 
state of the system, so that number of causes of visible 
randomness always highly exceeds number of its manifes¬ 
tations under time averaging even in most long realistic 
experiments. As the consequence, any particular exper¬ 
iment presents to researcher’s eyes its own unique as¬ 
sortment of relative frequencies (“probabilities”), or time 
rates, of random events composing a process under ob¬ 
servations. That is just the I/f-noise. 

Hence, being surprised at 1/f-noise is not more rea¬ 
sonable than being surprised at noise in general. The 
Nature needs 1/f-noise as expression of all inexhaustible 
resources of the Nature’s randomness in any particular 
“irreversible” processes as well as in originality of the 
wholly observed realization of our Universe’s evolution 
at all its time scales. A purely stochastic world, without 
1/f-noise, in which anything can be easily time-averaged, 
would be too tedious (and even, possibly, would repress 
a free will (H). 

Unfortunately, as we have seen above, 1/f-noise in¬ 
volves a “bad” statistics absolutely alien to the law of 
large numbers and resembling one what sometimes en¬ 
forces its observers, - for example, in [4l|, - to suspect 
action of mysterious “cosmic factors”. This fact signifi¬ 
cantly complicates theoretical tasks. 

Fortunately, although an influence from cosmos never 
is undoubtedly excluded, a source of randomness quite 
sufficient for 1 /f-noise creation is contained, - as we noted 
above, - already in so simple system as molecular Brow¬ 
nian particle interacting with ideal gas. And, gener¬ 
ally, - as we have demonstrated above, - a source of 
1/f-noise definitely exists in any medium which allows 
Brownian motion. Hence, one has every prospect of suc¬ 
cess in building and experimental verification of theory of 
1/f-noise and accompanying statistical anomalies start¬ 
ing from very usual Hamiltonians. 


We believe that the presented notes will induce some¬ 
body of interested readers to work in this intriguing area 


of statistical physics. 
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